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ABSTRACT 

An essential difference between two- 
dimensional and three-dimensional models 
of cycloidal rotors Is the presence of un- 
steady trailing cycloidal vorticies in the 
wake. The velocity Induced by these vor- 
ticies is the primary mechanism producing 
flow retardation for low span/radius 
ratio, finite blade number rotors. An 
Idealised rigid wake model of finite blade 
cycloidal rotors is used to Investigate 
some cycloidal rotor problems. 


INTRODUCTION 

During the past decade a rekindled 
Interest in energy conversion has stimula- 
ted the development of aerodynamic theories 
and numerical methods for analyzing rotors 
whose blades make cycloidal paths through 
the incident stream. Momentum theories 
and vortex theories appear in the litera- 
ture, In the case of vortex theories, 
some noteworthy models of the wake have 
been presented by Holme and by Wilson in 
References 1 and 2, 

Missing from these developments is 
Information about the velocity seen by the 
blades of a three-dimensional, finite 
blade number cycloidal rotor as It moves 
about Its orbit. Also missing is the 
answer to the question: ”What Is the 

optimum orbital load distribution?” There 
is good reason the latter question has 
gone unanswered: the problem Is complica- 

ted. In spite of research spanning seven 
decades, no theory or method having engin- 
eering usefulness has emerged for cycloidal 
propulsors. The power of present compu- 
ters may provide the means for pursuing 
this question further; it Is hoped that 
this work may be a step In that direction. 


DISCUSSION OF THE PROBLEM 

The load distribution on the blades 
of a cycloidal rotor has components due 

to: 

-flat plate angle of attack 
-camber 

-unsteady angle of attack (pitch) 
-unsteady motion normal to stream (heave) 
-circulation Induced by bound vorticity 
from other blades 

-circulation Induced from shed vorticity 
in the wake. 

Two-dimensional theory can be used to de- 
termine the load due to these components 
once the velocity at each orbital position 
and the character of the wake are determin- 
ed. 


An essential difference between the 
two-dimensional and three-dimensional 
models of cycloidal rotors is the presence 
of unsteady cycloidal vorticies in the 
wake. The velocity Induced by these vor- 
ticies is the primary mechanism producing 
flow retardation for finite and low span/ 
radius ratio, finite blade number rotors. 

It is for this reason that the author has 
chosen to investigate the velocity field 
of unsteady rigid cycloidal vorticies. 

The presence of unsteady bound and shed 
vorticity in the field between these vor- 
ticies Is recognized. The computation of 
the velocity induced by unsteady bound 
and shed vorticity is not within the scope 
of this paper, A comparison of the results 
of this paper with the results of the two- 
dimensional method of Mendenhall and Span- 
gler (Ref, 3) shows that the velocity in- 
duced by trailing cycloidal vorticies is 
several times that of the bound and shed 
vorticies. 

The wake xs assumed to be convected 
downstream at a constant velocity. It is 
further assumed that it does not deform. 
While such assumptions are not as defensi- 
ble for cycloidal rotors as they are for 
heiieoidal rotors, it can be argued that 
the Induced velocity seen by the blade is 
primarily Influenced by the portion of the 
cycloidal vortex filament which is nearest 
the blade. 

The validity of the model investigated 
may not prove useful for the complete 
range of tipspeeds and rotor loadings. 

One would expect it to be valid for very 
low tipspeed ratios where the model re- 
duces to a multielement airfoil system in 
rectilinear flow. One. would also expect 
it to be valid for lightly loaded rotors. 

It should be less satisfactory for a 
cycloidal propulsor whera a contracting 
wake will cause blades in the downstream 
half of the orbit to cut rolled vortex 
sheets originating In the upstream half of 
the orbit , The expanding wake of the 
turbine should minimize such effects. As 
with many other developments In the theory 
of fluid mechanics, the ultimate useful- 
ness can only be established when the com- 
plete induced velocity field is computed 
(including bound and shed vorticity effects) 
and when correlation with experimental 
data has been acheived. 


THE INDUCED VELOCITY 

We have in mind a discrete vortex 
lifting line model of a straight blade 
cycloidal turbine. We seek to evaluate 
the velocity induced at the blades of the 
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rotor by a system of cycloidal vorticies 
originating at different spanwise locations 
on the blades. One such vortex filament 
ant the coordinate system for this discus- 
sion are shown in Figure 1. The cycloidal 
curves followed by vorticies trailing from 
the i-th blade under rigid wake assumptions 
are parameterised in time by: 

x^ = R cos (Q-/^+f*>t ) + ut 

= R sin(9-f i -H^t) 
z i = z 

where the scalar quantity R is the radius 
of the blade orbit; u is a characteristic 
velocity through the rotor; 

*!= (i-l)2 y/N 

is the angular position of the i-th blade 
of an N-bladd rotor when blade 1 is at 
position 0 = 0;oj is the rotor angular 
velocity; and t is time (less than zero) 
representing conditions which orinigated 
in the past. 

According to the Biot-Savart law the 
velocity induced by a vortex filament is 
given by the following integral: 

V = 1 /\n f*p R x dl 
J I R| 3 

where P is A the circulation along the vortex 
filament; R is the vector from the point 
(x ,y , z ) where the velocity is computed 
to°th8 pjjint on the vortex fil- 

ament; dl is the unit tangent vector along 
the vortex filament directed away from the 
bound vortex. 

The circulation is a function of the 
spanwise position z, orbital position 0, 
and blade number i: 

Fi =A(z)B n cos(n9 + 0 n ) 

where Bp and 0 n give the variation in 
circulation around the orbit and A(z) is 
a spanwise load function. 

The velocity induced by this system 
has streamwise and sidewise components 
given by the following integrals: 

^6)=^- 'A(z)W4tf r 


APPLICATION OF THE MODEL 

A straight bladed cycloidal rotor 
having blade span/radius of l|/3 was div- 
ided into ten equal segments. A trailing 
cycloidal vortex system was modeled by 
eleven cycloidal vortex filaments origina- 
ting at the ends of these segments. The 
velocity induced at the midpoints of these 
segments was computed for a local tipspeed 
ratio Eol/u = 2.8£. The circulation on the 
segments was taken as the integrated aver- 
age over the length of the segment for an 
elliptical spanwise load distribution. 

The orbital load distribution was cosinu- 
soidal (Bn, 0 n = 0 for n>l). 

The velocity induced at the midpoints 
of the two segments nearest the midspan 
of a one blade cycloidal rotor 
Is shown In Figure 2. For the middle 60$ 
os the span the induced velocity is nearly 
constant. Near the tips it is not for two 
reasons. First, the discrete vortex lift- 
ing line model with equal segments does ' 
not yield the correct Induced velocity 
near the tips. In steady, rectilinear flow 
models the spacing Is usually modified near 
the tips. Secondly, there is no apriori 
guarantee that the circulation specified 
will result in the same pitch for the 
cycloidal vorticies that was used In com- 
puting the induced velocities. In fact, 
one must modify the pitch of the cycloidal 
vorticies and the repeat the calculation. 
Such Iteration is also required in solving 
the helicoidal rotor problem. The second 
Iteration has not been caried out in the 
work reported here. 

Another question arising from these 
results is: ”How do we determine the char- 
acteristic velocity at which the wake 
moves downstream?” We want, of course, to 
make a choice which gives the most accurate 
results. This question remains open, but 
In this work the velocity at 9 = 180 6 has 
been used. This is done because we will 
eventually want to use more terms In the 
cosine series for orbital circulation in 
the attempt to produce a more uniform 
induced velocity in the wake and the a 
more optimal rotor. 


cos (9-Yj_4wt+$) (z-z 0 )cos (0-y^-to>t) dt 
♦'-a, [(x-x 0 ) 2 + (y-y 0 ) z + (z-z o )*p/2 


v*j[9) = £~A(z)Buj/lyr f cos (9-^+1^+^) (z-Zg>) (sin(0-^+tJt ) -u/Ra>)dt 

J_ w {(x-x 0 ) 2 + (y-y 0 ) 2 + (z-z^P* 

For a straight blade rotor parallel to the 
z-axis we are not interested In the span- 
wise velocity component. These Integrals 
have been evaluated by numerical Integra- 
tion. 
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The results in Figure 2 form the 
basis on which the strength of the circu- 
lation and/or the incident stream may be 
varied to produce the local velocity at 
the rotor (u(l 80°)5 which matches the 
assumed local tipspeed ratio* If both 
are varied, a range of power coefficients 
is obtained. The largest then would be 
the optimum power coefficient for that 
local tipspeed ratio, or pitch, for that 
orbital and spanwise circulation distribu- 
%j ion , 

Figure 3 shows: 1 and 3 blade rotor 
results of varying the circulation and 
adjusting the inflow velocity until a 
local tipspeed ratio of 2.85 was acheived. 
These power coefficients are not quantita- 
tively exact since they do not account 
for induced velocities from shed and 
bound” vorticity, and are not derived from 
a completed iteration of spanwise loading 
and wake pitch. They do exhibit qualita- 
tive results showing an increase in obta- 
inable power coefficient with increased 
blade number and show that there is an 
optimal loading for each blade number, 

CONCLUDING REMARKS 

This method would be enhanced if 
a fast means of evaluating the induced 
velocity integral were found. Asymptotic 
series and a regression analysis of com- 
puted values seem reasonable approaches. 

The induced velocity components 
for hound vorticity and shed vorticity 
must be included in the analysis. 

The performance computation requires 
modifying the pitch of the trailing 
vorticles so that the wake mo tio n is con- 
sist ant with the prescribed inflow condi- 
tions and the induced velocities at the 
appropriate location. An iteration is 
required. 


The results of the inviscid theory 
must be reckonned with the reality of the 
viscous fluid effects. Correlation with 
experimental data is the ultimate test of 
a fluid dynamic theory or method. 

The method should be exercised in An 
attempt at the solution of the optimal 
orbital load distribution problem by in- 
cluding further terms in the circulation 
series. Assuming we can find the optimal 
load distribution, we must then see if it 
can be obtained by a technically feasible 
lifting surface. 

With minor extensions, this method 
could be applied to curved blade rotors 
of the Darrieus type. This is appealing 
because of the greater number of units 
which exist and the abundant theoretical 
and experimental date which is available 
for comparison. 
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